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ABSTRACT 

| We model an apparent gravitational lens system HDFS 2232509-603243in the Hub- 

Q\ [ ble Deep Field South. The system consists of a blue V = 25 mag arc separated by C//9 

from a red V = 22 mag elliptical galaxy. A mass distribution which follows the observed 
light distribution with a constant mass-to-light ratio can fit the arc component posi- 
tions if external shear is added. A good fit is also obtained with simple parameterized 
models, and all the models predict a fourth image fainter than the detection limit. The 
£Sj . inferred mass-to-light ratio is ~ 15 in solar units if the lens is at z = 0.6. The lens 

models predict a velocity dispersion of ~ 280 km s _1 which could be confirmed with 
t^. ' spectroscopy. 
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52 ■ 1. Introduction 

Detailed studies of gravitational lenses have provided a wealth of information on galaxies, both 
through modeling of individual lens systems (e.g., Schneider, Ehlers, & Falco 1992; Blandford &: 
Narayan 1992) and from the statistical properties of multiply imaged sources (e.g., Turner, Ostriker, 
& Gott 1984; Kochanek 1996). The extraordinary resolution of the Hubble Space Telescope (HST) 
has allowed for the first time the detection of lenses where both the source and the deflector are 
"normal" optically-selected galaxies. Ratnatunga et al. (1995) discovered two four-image "Einstein 
cross" gravitational lenses in the Groth-Westphal strip. Their mass models required ellipticities 
substantially above those of the light distributions, suggesting (as is commonly suspected in radio 
lenses) that some external shear is important to the lensing. 

The Hubble Deep Field (HDF, Williams et al. 1996), the deepest set of exposures taken with the 
HST, has enabled an unprecedented study of galaxy morphological evolution (e.g., Ellis 1998) and 
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the global star formation history (e.g., Madau, Pozzeti, & Dickinson 1998). Candidate gravitational 
lenses were reported by Hogg et al. (1996), who noted on the basis of radio surveys that several 
lenses are expected in the HDF (cf. Cooray, Quashnock, & Miller 1998). However, after further 
observations not one strong case for multiple imaging has been identified (Zepf, Moustakas, & Davis 
1997; Cohen 1998). 

The Southern Hubble Deep Field (HDF-S, Williams et al. 1998) is a second set of ultra-deep 
HST exposures, and several lenses are again expected. The best candidate lens, noted by the HDF- 
S team (Williams 1998), consists of a blue arc about CK'9 from an elliptical galaxy. In this Letter 
we report on a detailed study of this arc and argue that it is a bona fide multiply imaged source. 

2. Observations 

The WFPC2 HDF-S is a ~ 5 arcmin 2 field centered on 22 32 56.22 - 60 33 02.69 (J2000), 
imaged with the WFPC2 instrument on the HST in the F300W, FA50W, F606W, and F814H^ 
bandpasses for total exposure times of 1.4 x 10 5 , 1.0 x 10 5 , 8 x 10 4 , and 1.0 x 10 5 s respectively, 
reaching point-source sensitivities of roughly V = 29 mag (Williams et al. 1998). For the "version- 
1.0" reductions of the images, we adopt the Vega-relative calibrations supplied by the HDF-S team 
whereby 1 DN corresponds to F300W = 19.43 mag, FA50W = 22.02 mag, F606W = 22.90 mag, 
or 2^8141^ = 21.66 mag. The images clearly show a bright candidate lens system at pixel location 
(2982.7,2837.6). With the pixel scale of 0"0398 this is also 22 32 50.90 - 60 32 43.0 (J2000). This 
candidate lens system will be called "HDFS 2232509-603243" hereafter. 

The morphology of HDFS 2232509~603243is that of a V ~ 22 mag elliptical galaxy (hereafter 
"the elliptical") with a much bluer arc (hereafter "the arc") located / .'9 to the NW from the center 
of the elliptical, as shown in Figure |l|. There is also a faint point-like source (hereafter "the dot" ) 
located 0796 to the NE from the center of the elliptical. 

The elliptical is smooth and symmetrical, so images created by rotating its images through 
180° (centered on the pixel nearest the center of the elliptical) and subtracting them from their 
unrotated originals show essentially no sign of the elliptical. Photometry of the arc was performed 
in these "rotated-subtracted" images. The total F606W flux (given in Table [l]) was measured 
through a partial annulus of inner radius / .'517 and outer radius 1'.'194, extending over position 
angles 0° (north through east) through —108°. Colors were measured through a partial annulus 
of inner and outer radii / . / 677 and 1'.'035, over —18° through —90°. The flux and colors of the 
dot were measured through a / .'159 radius circular aperture. Because nearby pixels are correlated, 
signal-to-noise ratios for photometric measurements were found by re-binning the image into pixels 
4x4 times larger. 

The arc separates into four distinct components, labeled A, B, C and D, going clockwise 
around the arc. All have very similar colors. Components A, B, and D appear to have point-like 
centers with fuzz, while component C is more extended. Gaussians with the same FWHM as the 
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PSF (4 pixels) were fit to components A, B and D, and a gaussian of variable width was fit to 
component C, all in the FAhOW rotated-subtracted image. The best-fit positions and fluxes of 
these gaussians are given in Table [2|. 

In order to make arc-free images of the elliptical, the region of each image in the partial 
annulus centered on the elliptical with inner and outer radii of ff!517 and l'/194 and position angles 
0° through —108° was replaced with the corresponding region diametrically opposite (i.e., at a 
position angle rotated by 180°). These "arc-removed" images show virtually no evidence of the 
arc. The total F814W flux of the elliptical was measured out to its 2-cr (per pixel) isophote. The 
colors were measured through 0'.'915 radius circular apertures, and are given in Table [l]. 

A de Vaucouleurs profile is a very good fit to the azimuthally averaged radial profile of the 
elliptical in the "arc-removed" images. The effective radius is r e = (X'31 and the F814W-band 
surface brightness at the effective radius is 21.0 mag in each arcsec 2 . The isophote ellipticity 
orientation varies with radius. The outer isophotes (radii > Of! 8) of the elliptical have ellipticity 
~ 0.2 and are elongated along +40°. The inner isophotes (radii < 0f.'3) also have ellipticity ~ 0.2 
but are elongated along —40°, i.e., almost orthogonal. Isophotes intermediate to these two regimes 
are fairly circular. 

Comparison of the colors of the elliptical with HDF galaxies with known redshifts suggests that 
its redshift is roughly z = 0.6 (cf. the photometric redshift of 0.5 of Gwyn (1998)). This suggests 
a total luminosity^] of 8.7 x 1O 9 /i~ 2 L0 in the rest-frame B-band, or roughly L*. The luminosity is 
a factor of 3 larger if z = 0.8, and a factor of 6.5 smaller if z = 0.3. 

3. Lens models 

We consider several models for the lensing mass distribution, simple parameterized models as 
well as a constant mass-to-light ratio model. The image positions are not very sensitive to the 
radial profile of the lens, so for the parameterized models we take it to be singular and isothermal. 
Asymmetry is required and we consider two types, external shear or an elliptical galaxy. We assume 
that the arc components A, B and D are three images of a common source, with A and B a merging 
pair. Component C can easily be produced if the source is extended or multi-component. We do 
not include as constraints component C or the observed flux ratios among components A, B and D, 
since these are sensitive to the detailed flux distribution of the source. We discuss the dot below. 

We consider a lens at redshift zl and a source at zs- Then the lens mass distribution determines 
a convergence k and a potential ip (e.g., Schneider et al. 1992; Blandford & Narayan 1992). The 
singular isothermal sphere (SIS) is defined by k = b/[2r]. We add to the potential the term 
ipry = —^"fr 2 cos 2(9 — 9y) , which defines an external shear of magnitude 7 and direction # 7 . If the 
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shear is due to an axisymmetric galaxy then it is located at an angle # 7 from the lens galaxy. We 
denote by SIS+7 the SIS model with external shear. 

A singular isothermal elliptical mass distribution (SIEMD) with an axis ratio a and a major 
axis along the y-axis is given by k = b/[2y/ x 2 + y 2 a 2 ]. More generally we rotate the major axis by 
an angle 6 e . The lensing properties of the SIEMD have been studied by Kassiola & Kovner (1993) 
and by Kormann, Schneider, & Bartelmann (1994). 

We fit these two lens models, the SIS+7 model and the SIEMD, to the three observed image 
positions, with results shown in Table [| Both lens models have an ellongated diamond caustic 
with four cusps. We indicate a source near a cusp on the major axis by [maj] and one near a minor 
cusp by [min]. Since the observed position B is farther from the center of the lens than A or D, 
the arc is more naturally produced by a minor cusp for which only the middle image of the three 
is outside the critical curve. 

Interestingly, the orientation of the SIEMD [min] model agrees with the observed orientation of 
the elliptical's outer isophotes, while the orientation of the SIEMD [maj] model agrees with that of 
the inner isophotes. However, among the simple parameterized models only the SIS+7[min] model 
fits the data well (see Table || and Figure ||). This model requires a substantial shear (7 = 0.26) 
produced by objects lying in the NE or SW directions. To compute the mass-to-light ratio for 
various models, we fix a circle about the galaxy center of radius 07915. The luminosity of the 
galaxy within this circle is 7.6 x 10 9 /i _2 Lq if zl = 0.6. If, e.g., z$ = 1.5, then the mass within 
the same circle of the SIS is M = 2.1 x 10 11 /i _1 (6/1")Mq. This yields a mass-to-light ratio of 
27 h (6/1"), and a velocity dispersion a = 270-i/fr/l" km s _1 . 

As a second approach to modeling, we utilize the resolved image of the lens galaxy (smoothed 
with a 4-pixel FWHM gaussian) to construct a model with a constant mass-to-light ratio (denoted 
M/L). The pure M/L model fails but together with external shear it produces an excellent fit, 
albeit with a large shear. Even a small shear such as 0.05 produces a substantially better fit than 
with no shear (see Table |3| and Figure ||). All the M/L models produce a minor cusp configuration. 
The parameter R for the M/L models gives the mean k in a circle of radius 07915 about the galaxy 
center. If zl = 0.6 and zs = 1.5 then the mass within this circle is M = 1.9 x lO 11 ^ 1 RMq, 
which yields a mass-to- light ratio of 25 hR. If we assume a small shear then R ~ 1 (Table ||), and 
the mass-to-light ratio is on the high side of the range of other lens galaxies (Keeton, Kochanek, & 
Falco 1998). Note, however, that we only have an estimated photometric redshift for the lens and 
a guess for the source redshift. The mass-to-light ratio changes to 14 hR for {zl = 0.8, zs = 2.0) 
and to 57 hR for {zl = 0.3, zs = 1.0). 

No fourth image is detected down to the detection threshold of about 28.9 mag in F606W. 
This is consistent with the models, which predict a fourth image fainter than the third by a factor 
of 3 or 4. The M/L models are non-singular and predict a fifth image fainter by an additional factor 
of ~ 100. The best-fitting models also locate the image of the cusp near the observed position of 
component C, so this component can be produced if the source is extended. 
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Thus far we have not included the dot, assuming it to be a separate background source. We 
consider also the possibility that the dot is the third image corresponding to components A and B. 
The types of lens models used above can fit these alternative image positions with similar x 2 values 
as in Table 3, but they fail since they predict an observable arc which includes counter images of 
components C and D and extends South-East from the dot. More elaborate lens models could work 
but only if the sources of components C and D are small (subpixel) and lie right on the caustic. 

4. Discussion 

We have modeled a candidate gravitational lens, HDFS 2232509-603243, found in the HDF-S. 
The image positions can be fit in a minor cusp configuration by simple lens models which include 
cllipticity or shear, or with shear added to a model with a constant mass-to-light ratio. However, 
an exact fit is obtained only with a substantial shear produced by objects lying to the NE or SW. 
No obvious candidate galaxy lies within 5" of the lens, but there are several galaxies ~ 20" to the 
SW and others ~ 35" to the NE with photometric redshifts in the range 0.5-0.7 (Gwyn 1998). 
Spectroscopic redshifts will determine the possible presence of galaxy groups. 

The lensing galaxy is expected to have a redshift z ~ 0.6, which gives it a luminosity of ~ L*. 
Lens models imply a mass-to-light ratio of ~ 25 h, and a velocity dispersion of ~ 280 km s -1 . 
If the elliptical falls on the fundamental plane, a central line-of-sight velocity dispersion of only 
~ 100 km s~ x is expected if z = 0.6 (Jorgensen, Franx & Kjaergaard 1996), which suggests that 
z ~ 0.8 is more likely. 

Showing that the arc and the lens galaxy are at different redshifts would rule out the possibility 
that the arc is a tidally disrupted satelite. On 10-m-class telescopes, it is straightforward to get 
spectroscopic redshifts of galaxies at 22 mag, even elliptical galaxies which tend to have absorption- 
feature dominated spectra (Cohen et al. 1998). Some elliptical galaxies at this flux level even have 
measured central velocity dispersions (van Dokkum et al. 1998; Pahre 1998). Thus the prospects 
for follow-up spectroscopy on the elliptical are good. It is common to obtain redshifts for sources 
at 25 mag, even some which are close to bright nearby sources (Steidel et al. 1996). A redshift for 
the arc is likely obtainable if it is at z < 1.3 since its blue color suggests that it will show strong 
emission lines in a visual spectrum (Hogg et al. 1998). The arc is unlikely to be at redshift z > 2.3 
since it shows significant F300W flux. 

It is a pleasure to thank Bob Williams and the entire HDF-S team for taking, reducing, and 
making public the beautiful images of the HDF-S. Scott Buries, Wayne Hu and Inger Joergensen 
provided very useful help in short order. We thank the referee Emilio Falco for useful comments. 
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Fig. 1. — A region taken from the public release of the HDF-S (Williams et al. 1998). It shows the 
lens galaxy, the arc to its NW, and the dot to its NE. North is up and East is to the left (except 
for a 0.5° rotation). The figure is centered on the lens galaxy and measures 3'.'2 on a side. 
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Fig. 2. — Image and source positions for the SIS+7[min] model (column 3 of Table ||). Four 
image positions (o) are indicated along with the source position (•) and the observed positions of 
arc components A, B and D (*). Also shown are the tangential caustic (solid) and critical curve 
(dotted). The figure is centered on the lens galaxy, with the a-direction west and the (^-direction 
north. 
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Fig. 3. — Same as Figure 2, but for the M/L+7 model with 7 = 0.05 (last column of Table ||[) 
There are five images, although the one near the lens center is extremely faint. 
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Table 1. Photometry 



component 


flux a 




F300W 


- FA50W 


F450W 


- F606W 


F606W 


- F8UW 




(mag) 


[snr] b 


(mag) 


[snrl,snr2] c 


(mag) 


[snrl,snr2] 


(mag) 


[snrl,snr2] 


elliptical 


F8UW = 20.52 


[950] 


0.31 


[9.2,99] 


1.86 


[99,610] 


1.53 


[610,1200] 


arc 


^6061^ = 25.19 


[41] 


-0.96 


[11,37] 


0.24 


[37,51] 


0.40 


[51,37] 


dot 


F606W = 28.70 


[6.2] 




[-1.5,3.9] 


0.39 


[3.9,6.2] 


-0.05 


[6.2,2.9] 



a Different components have flux measurements in different bands and through different focal- 
plane apertures; see text for details. 

b Signal-to-noise ratios (snr) measure the flux through the aperture relative to the expected sky 
noise rms through that same aperture. See text for aperture definitions. 

c The snr values snrl and snr2 are for the two bandpasses involved in the color measurement, 
bluer first. 



Table 2. Astrometry a and Photometry 



sub-component 


Aa 


A5 


FA50W h 




(arcsec) 


(arcsec) 


(mag) 


elliptical 










arc-A 


+0.406 


+0.777 


27.55 


arc-B 


+0.568 


+0.721 


27.43 


arc-C 


+0.711 


+0.478 


27.08 


arc-D 


+0.824 


+0.192 


28.10 


dot 


-0.778 


+0.567 





Positions are given relative to 22 32 50.9017 -60 32 43.009 (J2000), with the a-direction west 
and the 5-direction north. We estimate position errors of / . / 02 from fitting an artificial image. 

b The magnitudes of the arc components are based on the fluxes of the fitted gaussians. The 
fluxes do not add up to the total arc flux because the magnitudes exclude some diffuse emission. 
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Table 3. Best fit model results 



Parameter 


SIEMD[min] 


SIEMD[maj] 


SIS+7[min] 


SIS+7[maj] 


M/L 


M/L+7 


M/L+- 


b or R 


0785 


0"84 


i'.'os 


0"87 


1.10 


1.40 


0.998 


a or 7 


0.031;^ 


0.961;^ 


o.26i;8* 


0.005i;^ 5 




0.52i;^ 


0.05 


9 e or 9j 


41?2±1;| 


-48?6i|? 


42? 5t- 7 7 


-49?0iy 




43?5±;| 


53? 4 


#dof 


1 


1 


1 


1 


3 


1 


1 


t 


5.1 


5.6 


0.26 


5.8 


97 


0.40 


4.5 



Note. - The best fits are formally at a — > for the SIEMD[min] model, at a — ► 1 for the 
SIEMDfmaj] model, and at 7 — > for the SIS+7[maj] model. The overall best fit for the M/L+7 
model is given in column 6, while column 7 gives the best fit for this model with the restriction 
7 = 0.05. Where given, parameter ranges indicate la uncertainties defined by Ax 2 = X 2 - All 
angles are measured north through east. 



